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Abstract 

Introduction: In recent genome-wide association studies for psoriatic arthritis (PsA) and psoriasis vulgaris, common 
coding variants in the TRAF3IP2 gene were identified to contribute to susceptibility to both disease entities. The risk 
allele of p.AsplOAsn (rs33980500) proved to be most significantly associated and to encode a mutant protein with 
an almost completely disrupted binding property to TRAF6, supporting its impact as a main disease-causing variant 
and modulator of IL-17 signaling. 

Methods: To identify further variants, exons 2-4 encoding both known TNF-receptor-associated factor (TRAF) 
binding domains were sequenced in 871 PsA patients. Seven missense variants and one three-base-pair insertion 
were identified in 0.06% to 1.02% of alleles. Five of these variants were also present in 931 control individuals at 
comparable frequency. Constructs containing full-length wild-type or mutant TRAF3IP2 were generated and used to 
analyze functionally all variants for TRAF6-binding in a mammalian two-hybrid assay. 

Results: None of the newly found alleles, though, encoded proteins with different binding properties to TRAF6, or 
to the cytoplasmic tail of the IL-17-receptor a-chain, suggesting that they do not contribute to susceptibility. 

Conclusions: Thus, the 7/?/\F3/P2-variant p.AsplOAsn is the only susceptibility allele with functional impact on 
TRAF6 binding, at least in the German population. 



Introduction 

Recent genome-wide association studies (GWASs) for 
psoriatic arthritis and psoriasis vulgaris have been quite 
successful and revealed more than 25 new susceptibility 
loci [1-7] characterizing psoriatic arthritis (PsA), mainly 
as an immune-mediated disease that is also affected by 
barrier proteins. One new PsA and psoriasis-susceptibil- 
ity gene, TRAF3IP2, codes for the adaptor protein ACT1 
(nuclear factor-^B activator 1), a regulator of the NF-^B 
pathway involved in IL-17 signaling. We were able to 
identify the common coding variant p.AsplOAspn 
(rs33980500) as the disease-causing allele with evidence 
from association findings at single nucleotide poly- 
morphism (SNP) and haplotype levels, as well as binding 
studies of ACT1 with its interaction partner TRAF6 [2] 
(Figure 1). However, despite the exponentially increased 
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knowledge of PsA and psoriasis susceptibility, the sum 
of the currently identified genetic risk factors does not 
explain the entire genetic contribution to disease. 

At the common susceptibility gene IL23R, coding for 
interleukin-23 receptor for Crohn disease, low-frequency 
variants have been identified as protective factors, and 
their effects turned out to be larger than those of common 
variants and accounted for a greater portion of genetic 
contribution [8]. To find out whether low- frequency cod- 
ing variants in TRAF3IP2 contribute to PsA susceptibility, 
we sequenced exons encoding TRAF-binding domains 
(Figure 1) [9] and studied their association with PsA as 
well as the functional consequences of the respective 
amino acid changes on interaction with TRAF6 and the 
interleukin-17 receptor a (IL17-RA). 

Materials and methods 

Study population 

We analyzed a case-control study comprising 871 PsA 
patients. A subset of 764 patients was described previously 
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Figure 1 Interleukin (IL)-1 7A/IL1 7F binds to the heterodimeric 
IL-17R and leads to its signaling. Signaling includes recruitment 
of nuclear factor-^B activator 1 (ACT1) via SEFIR domains (red ovals). 
Activation of ACT 1 encoded by TRAF3IP2 allows binding to TRAF6 
via TNF receptor-associated factor (TRAF)-binding domains (yellow 
ovals); this activates NF-^B as well as proinflammatory cytokines. 



[2]. All patients fulfilled the recently defined CASPAR (for 
ClASsification of Psoriatic ARthritis) criteria [10] and were 
recruited between 2002 and 2011 exclusively by board-cer- 
tified rheumatologists at five different rheumatology 
departments in Germany. In comparison to the previously 
described cohort, clinical characteristics were very similar: 
the mean (SD) age of onset for PsV was 28.3 ± 13.3 years; 
60.8% of the patients were men. For 78% of the patients, 
the diagnosis of PsA was made > 3 years before recruit- 
ment, and 96% of patients had a skin involvement > 3 
years before study inclusion. Peripheral joint involvement 
was detectable in the majority of cases (827 or 94.8%); this 
was oligoarticular in 183 patients and polyarticular in 644 
(21.0% and 73.8% of the entire cohort, respectively). Diag- 
nosis of spinal involvement was based on symptoms of 
inflammatory back pain, characteristic clinical signs of 
restricted vertebral movement and/or sacroiliac pain on 
physical examination, and a subsequent confirmation by 
radiographic signs of sacroiliitis or spondylitis or both. 
Spinal involvement was observed in 204 patients, account- 
ing for 23.4% of the PsA cohort. In these patients, sacroilii- 
tis or spondylitis or both were partly associated with 
concomitant peripheral joint disease. 



The 931 German control individuals were healthy 
blood donors and the same as previously reported [2]. 

Ethical approvals for all German studies were obtained 
from the Ethics Committees of the medical faculties of the 
Universities of Erlangen, Munster, and Frankfurt/Main. 
All probands gave written informed consent. The investi- 
gations were conducted according to the Declaration of 
Helsinki principles. 

Sequencing 

All PsA patients were sequenced for exons 2-4 coding 
for TRAF-binding domains. Mainly intron-based pri- 
mers for four amplicons were designed with Primer 3 
(v.0.4.0) [11]. In general, Sanger sequencing was per- 
formed, as described before [12], with the following 
changes: thermocyclers were Thermo Multiblock (Ther- 
moFisher Scientific, Ulm, Germany); the robotic system 
for PCR (polymerase chain reaction) clean-up as well as 
for sequencing reactions was Beckman-Coulter Biomek 
Nxp (kit for PCR clean-up: Agencourt AMPure; kit for 
sequencing reactions: CleanSEQ; Beckman-Coulter, Kre- 
feld, Germany). Sequences were analyzed with the Soft- 
ware Sequencher versus 4.10.1 (Gene Codes, Ann 
Arbor, MI, USA). Genotyping rates of the four ampli- 
cons were between 99.8% and 100%. Variants were 
named according to the usual naming conventions with 
regard to TRAF3IP2 reference sequence NM_147686 
(coding for 565 amino acids). 

Analysis of variants in controls 

Control individuals were genotyped for all missense 
variants as well as the 3bp-insertion with self-designed 
TaqMan assays (Life Technologies, Carlsbad, CA, USA), 
as previously described [2]. Missing genotypes or carriers 
of low-frequency variants were resequenced, resulting in a 
genotyping rate of 100% for seven variants and 99.8% for 
the further one. 

Statistics 

Fisher Exact test, as implemented in R (version 2.13.1) 
[13], was used for comparison of frequency distributions 
between cases and controls. 

Protein-binding analyses in mammalian two-hybrid 
system 

A mammalian two-hybrid assay (Stratagene) was applied 
to analyze the interaction of wild-type and mutant protein 
ACT1 of TRAF3IP2 with TRAF6 and IL-17RA, as 
described before [2]. As a change to previous analyses, 
full-length TRAF3IP2 mRNA (coding for amino acid resi- 
dues 1 through 565) was amplified from reverse-tran- 
scribed cDNA of the chondrocyte cell line T/C28a4 by 
using S'-TGAATTCATGAACCGAAGCATTCCTGTG-S' 
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and 5^TGCGGCCGCTCACAAGGGAACCACCTGA 
AG-3\ The cytoplasmic tail of IL-17RA (amino acid resi- 
dues 343-667, NM_014339) was amplified from reverse- 
transcribed cDNA of the human macrophage cell line 
THP-1 by using S'-AGGATCCATGACCTGGAGGC- 
TAGCTGGGCCTGGAA-3' and S'-AGCGGCCGCT- 
CAGGTGTGGAGGGGCTGCGGCGCTGGCTGA-3', 
and cloned into the BamHl/Notl site of the pCMV-BD 
bait vector. HEK 293 cells (1 x 10 4 ) were grown in 96- well 
white tissue-culture plates (Greiner) for 24 hours, and the 
bait and prey vectors (25 ng each) were co-transfected 
with a firefly luciferase reporter (250 ng) by using jetPEI 
reagents according to the suppliers instructions (PeqLab). 
Transfection efficiency was controlled by co-transfecting a 
renilla luciferase plasmid (5 ng, pRL-TK vector; Promega). 
At 48 hours after transfection, luciferase activity of quad- 
ruplicate wells was measured by using the Dual Glo Assay 
System (Promega) on a Mithras LB940 plate reader (Bert- 
hold Technologies). Assays were performed > 5 times. 

For functional testing of the two common as well as the 
newly identified rare coding variants of TRAF3IP2 [2], 
constructs were generated by site-directed mutagenesis by 
using the QuickChange site-directed mutagenesis kit 
(Stratagene). Quality control of all vectors was performed 
with DNA sequencing. The firefly luciferase values were 
normalized to the Renilla values. Statistical significance 
was determined by using an unpaired Student t test. 

Results 

Sanger sequencing in 871 PsA patients revealed eight low- 
frequency variants: one 3-bp insertion and seven missense 



variants, one of them being the known SNP rs61756667 
(minor allele frequency of 0.10% in the database), as indi- 
cated in Table 1. Allele frequencies ranged between 0.06% 
and 1.03% in PsA cases and 0 to 0.54% in control indivi- 
duals, respectively. Distributions of allele frequency for 
single low-frequency variants as well as the sum of the 
eight variants (P = 0.493) did not show significant differ- 
ences between cases and controls. 

Because our case-control study of -900 individuals was 
of limited power to exclude association to very-low- 
frequency variants (for example, found on one allele of 
one patient), we investigated the impact of amino acid 
changes encoded by the new TRAF3IP2-vaiciants on func- 
tional properties of the ACT1 protein also by mammalian 
two-hybrid technology. None of the amino acid muta- 
tions in the full-length constructs of the new TRAF3IP2 
coding variants exhibited any influence on TRAF6 bind- 
ing (Figure 2). Also, TRAF6 interaction of the common 
coding variant rsl3 190932 (p.Arg74Trp) was comparable 
to wild-type, in agreement with earlier results obtained in 
assays using shortened TRAF3IP2 constructs [2] . In con- 
trast, the analysis of a full-length TRAF3IP2 construct, 
containing the earlier described disease-causing variant 
rs33980500 (p.AsplOAsn), revealed a nearly completely 
disrupted interaction with TRAF6 (-90%), thereby con- 
firming previous analyses. Moreover, the ability to inter- 
act with the cytoplasmic IL-17RA domain did not differ 
between ACT1 proteins encoded by either wild-type 
TRAF3IP2 or by any of its low-frequency variants. In this 
respect, p.AsplOAsn did not cause different results com- 
pared with other TRAF3IP2 variants. 



Table 1 Association analysis of new variants 



Coding position 


Protein position 


Position hg19 


Allele 


871 PsA patients 


931 control individuals 


P (Fisher Exact) 










n 


% 


n 


% 




C.105J06 InsACC 


p.Pro35dup 


111913185 


Ins 


2 


0.11 


3 


0.16 


1 








WT 


1,740 


99.89 


1,859 


99.84 




C.281G > A (rs61 756667) 


p.Ser94Asn 


111913009 


A 


1 


0.06 


4 


0.21 


0.376 








G 


1,737 


99.94 


1,858 


99.79 




C.350C > T 


p.Alal 17Val 


111912940 


T 


1 


0.06 


0 


0 


0.483 








C 


1,737 


99.94 


1,858 


100.00 




C.649C > A 


p.Pro217Thr 


111912641 


A 


18 


1.03 


10 


0.54 


0.127 








C 


1,724 


98.97 


1,852 


99.46 




C.691C > T 


p!eu231Phe 


111912599 


T 


2 


0.11 


1 


0.05 


0.613 








C 


1,740 


99.89 


1,861 


99.95 




C.746A > G 


p.Gln249Arg 


111912544 


G 


1 


0.06 


0 


0 


0.483 








A 


1,741 


99.94 


1,862 


100.00 




C.1058G > C 


p.Gly353Ala 


111896989 


C 


1 


0.06 


0 


0 


0.483 








G 


1,741 


99.94 


1,862 


100.00 




C.1184A > G 


p.Asn395Ser 


111896863 


G 


3 


0.17 


7 


0.38 


0.346 








A 


1,739 


99.83 


1,855 


99.62 





WT, wild type. Names, position on chromosome 6, and count/frequency of wild-type alleles/low-frequency alleles of new variants in PsA patients and control 
individuals are shown, as well as the P value of Fisher Exact test. 
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Figure 2 Binding of ACT1 and its mutant variants to TNF receptor-associated factor (TRAF)6 and interleukin (IL)-17RA. Interaction of 
wild-type TRAF 3 1 P2/ ACT 1 and its mutant variants with either TRAF6, upper panel, or IL-17RA, lower panel, was determined by mammalian two- 
hybrid system. The TRAF3IP2 mutant p.AsplOAsn showed a significantly reduced binding to TRAF6. Common coding variants indicated by 
underlined and italicized letters. *P < 0.005, when compared with wild-type TRAF3IP2. 



Discussion 

Our in-depth sequencing of TRAF-binding domains in 
TRAF3IP2 revealed a number of new low-frequency var- 
iants in PsA patients occurring at a similar frequency in 
control individuals. This finding is not unexpected, 
because, for example, newly identified low-frequency var- 
iants in most of the susceptibility genes identified in 
GWAS and sequenced for Crohn disease [8] or type I 
diabetes [14] were not associated with disease. Further 
evidence that our newly identified coding variants do not 
cause disease was derived from functional analyses: none 
of the variants leading to amino acid changes in the 
respective ACT1 proteins affected TRAF6 interaction. 
Reassuringly, we could confirm previous results of a dra- 
matically reduced TRAF6 binding for the ACT1 variant 
harboring the disease-causing variant p.AsplOAsn under 
assay conditions analyzing full-length TRAF3IP2-tr2in- 
scripts (Figure 2). Binding of ACT1 to the cytoplasmic 
domain of IL-17RA is mediated via a C-terminal SEFIR 
domain that remained conserved in wild-type configuration 



in all investigated TRAF3IP2 constructs (Figures 1 and 2). 
Accordingly, a rather limited a priori likelihood existed to 
detect major differences in IL-17RA-binding properties of 
the investigated mutant TRAF3IP2 constructs, although an 
indirect modulatory effect of conformational changes 
induced by the distantly located mutations in the TRAF- 
binding domains could not be entirely excluded. However, 
our studies did not reveal an altered IL-17RA interaction 
for any of the nine investigated coding variants. 

Conclusions 

Our data indicate that the common variant rs33980500 
(p.AsplOAspn) is the only risk allele for PsA with func- 
tional impact on TRAF6 binding, at least in the German 
population. 

Abbreviations 

ACT1: nuclear factor-KB activator 1; CASPAR: ClASsification of Psoriatic 
Arthritis; GWAS: genome-wide association studies; IL-17RA: interleukin-1 7 
receptor A; NF-kB: nuclear factor-KB; PCR: polymerase chain reaction; PsA: 



Bohm et al. Arthritis Research & Therapy 2012, 14:R84 
http://arthritis-research.eom/content/14/2/R84 



Page 5 of 5 



psoriatic arthritis; PsV: psoriasis vulgaris; SNP: single-nucleotide 
polymorphism; TRAF: TNF receptor-associated factor; TRAF3: TNF receptor- 
associated factor 3; TRAF6: TNF receptor-associated factor 6; TRAF3IP2: 
TRAF3-interacting protein 2. 

Acknowledgements 

We are grateful to all patients and control probands for participation in this 
study. We thank Kristin Krause and Petra Badorf for excellent technical 
assistance. 

The work was supported in part by a grant from the Interdisciplinary Centre 
for Clinical Research (J1; to UH) of the University of Erlangen-Nuremberg, 
Germany, as well as the German Federal Ministry of Education and Research 
[Bundesministerium fur Bildung und Forschung (BMBF); project ANCYLOSS; 
to AR and ArthroMark (project 4, 01 EC 1009C); to HB]. 

Author details 

1 Divisions of Rheumatology, Department of Internal Medicine II, Johann 
Wolfgang Goethe University, Theodor-Stern-Kai 7, 60590 Frankfurt am Main, 
Germany, institute of Human Genetics, University of Erlangen-Nuremberg, 
Schwabachanlage 10, 91054 Erlangen, Germany. 

Authors' contributions 

HB, AR, and UH designed the study and worked out its concept. UH 
planned and performed the genotyping, and SU and UH analyzed genetic 
data. BB and HB performed functional studies. BB, HB, MA, FB, and UH 
recruited subjects and collected phenotypic data. UH, BB, and HB wrote a 
first draft of the manuscript; all authors reviewed and approved the 
manuscript. 

Competing interests 

The authors declare that they have no competing interests. 

Received: 2 November 2011 Revised: 22 March 2012 
Accepted: 18 April 2012 Published: 18 April 2012 

References 

1. Ellinghaus E, Ellinghaus D, Stuart PE, Nair RP, Debrus S, Raelson JV, 
Belouchi M, Fournier H, Reinhard C, Ding J, Li Y, Tejasvi T, Gudjonsson J, 
Stoll SW, Voorhees JJ, Lambert S, Weidinger S, Eberlein B, Kunz M, 
Rahman P, Gladman DD, Gieger C, Wichmann HE, Karlsen TH, Mayr G, 
Albrecht M, Kabelitz D, Mrowietz U, Abecasis GR, Elder JT, et al: Genome- 
wide association study identifies a psoriasis susceptibility locus at 
TRAF3IP2. Not Genet 2010, 42:991-995. 

2. Huffmeier U, Uebe S, Ekici AB, Bowes J, Giardina E, Korendowych E, 
Juneblad K, Apel M, McManus R, Ho P, Bruce IN, Ryan AW, Behrens F, 
Lascorz J, Bohm B, Traupe H, Lohmann J, Gieger C, Wichmann HE, 
Herold C, Steffens M, Klareskog L, Wienker TF, Fitzgerald O, Alenius GM, 
McHugh NJ, Novelli G, Burkhardt H, Barton A, Reis A: Common variants at 
TRAF3IP2 are associated with susceptibility to psoriatic arthritis and 
psoriasis. Not Genet 2010, 42:996-999. 

3. Nair RP, Duffin KC, Helms C, Ding J, Stuart PE, Goldgar D, Gudjonsson JE, 

Li Y, Tejasvi T, Feng BJ, Ruether A, Schreiber S, Weichenthal M, Gladman D, 
Rahman P, Schrodi SJ, Prahalad S, Guthery SL, Fischer J, Liao W, Kwok PY, 
Menter A, Lathrop GM, Wise CA, Begovich AB, Voorhees JJ, Elder JT, 
Krueger GG, Bowcock AM, Abecasis GR: Genome-wide scan reveals 
association of psoriasis with IL-23 and NF-kappaB pathways. Not Genet 
2009, 41:199-204. 

4. Zhang XJ, Huang W, Yang S, Sun LD, Zhang FY, Zhu QX, Zhang FR, 
Zhang C, Du WH, Pu XM, Li H, Xiao FL, Wang ZX, Cui Y, Hao F, Zheng J, 
Yang XQ, Cheng H, He CD, Liu XM, Xu LM, Zheng HF, Zhang SM, Zhang JZ, 
Wang HY, Cheng YL, Ji BH, Fang QY, Li YZ, Zhou FS, et al: Psoriasis 
genome-wide association study identifies susceptibility variants within 
LCE gene cluster at 1q21. Not Genet 2009, 41:205-210. 

5. Strange A, Capon F, Spencer CC, Knight J, Weale ME, Allen MH, Barton A, 
Band G, Bellenguez C, Bergboer JG, Blackwell JM, Bramon E, Bumpstead SJ, 
Casas JP, Cork MJ, Corvin A, Deloukas P, Dilthey A, Duncanson A, Edkins S, 
Estivill X, Fitzgerald O, Freeman C, Giardina E, Gray E, Hofer A, Huffmeier U, 
Hunt SE, Irvine AD, Jankowski J, et al: A genome-wide association study 
identifies new psoriasis susceptibility loci and an interaction between 
HLA-C and ERAP1. Not Genet 2010, 42:985-990. 



10. 



14. 



Stuart PE, Nair RP, Ellinghaus E, Ding J, Tejasvi T, Gudjonsson JE, Li Y, 
Weidinger S, Eberlein B, Gieger C, Wichmann HE, Kunz M, Ike R, Krueger GG, 
Bowcock AM, Mrowietz U, Lim HW, Voorhees JJ, Abecasis GR, 
Weichenthal M, Franke A, Rahman P, Gladman DD, Elder JT: Genome-wide 
association analysis identifies three psoriasis susceptibility loci. Nat Genet 
2010, 42:1000-1004. 

Sun LD, Cheng H, Wang ZX, Zhang AP, Wang PG, Xu JH, Zhu QX, Zhou HS, 
Ellinghaus E, Zhang FR, Pu XM, Yang XQ, Zhang JZ, Xu AE, Wu RN, Xu LM, 
Peng L, Helms CA, Ren YQ, Zhang C, Zhang SM, Nair RP, Wang HY, Lin GS, 
Stuart PE, Fan X, Chen G, Tejasvi T, Li P, Zhu J, et al: Association analyses 
identify six new psoriasis susceptibility loci in the Chinese population. 
Not Genet 2010, 42:1005-1009. 

Momozawa Y, Mni M, Nakamura K, Coppieters W, Aimer S, Amininejad L, 
Cleynen I, Colombel JF, de Rijk P, Dewit O, Finkel Y, Gassull MA, 
Goossens D, Laukens D, Lemann M, Libioulle C, O'Morain C, Reenaers C, 
Rutgeerts P, Tysk C, Zelenika D, Lathrop M, Del-Favero J, Hugot JP, de 
Vos M, Franchimont D, Vermeire S, Louis E, Georges M: Resequencing of 
positional candidates identifies low frequency IL23R coding variants 
protecting against inflammatory bowel disease. Nat Genet 201 1, 43:43-47. 
Hunter CA: Actl-ivating IL-17 inflammation. Nature Immunol 2007, 
8:232-234. 

Taylor W, Gladman D, Helliwell P, Marchesoni A, Mease P, Mielants H: 
Classification criteria for psoriatic arthritis: development of new criteria 
from a large international study. Arthritis Rheum 2006, 54:2665-2673. 
Primer 3 vs. 0.4.0, Pick primers from a DNA sequence, [http://frodo.wi.mit. 
edu/]. 

Huffmeier U, Lascorz J, Traupe H, Bohm B, Schurmeier-Horst F, Stander M, 
Kelsch R, Baumann C, Kuster W, Burkhardt H, Reis A: Systematic linkage 
disequilibrium analysis of SLC12A8 at PSORS5 confirms a role in 
susceptibility to psoriasis vulgaris. J Invest Dermatol 2005, 125:906-912. 
R = language and environment for statistical computing and graphics. 
[http://www.r-project.org/]. 

Nejentsev S, Walker N, Riches D, Egholm M, Todd JA: Rare variants of 
IFIH1, a gene implicated in antiviral responses, protect against type 1 
diabetes. Science 2009, 324:387-389. 



doi:10.1186/ar3807 

Cite this article as: Bohm et al.: Identification of low-frequency TRAF3IP2 
coding variants in psoriatic arthritis patients and functional 
characterization. Arthritis Research & Therapy 2012 14:R84. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at 
www.biomedcentral.com/submit 



o 



BioMed Central 



